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ABSTRACT 

Establishing relations between global stellar parameters and asteroseismic quantities 
can help improve our understanding of stellar astrophysics and facilitate the interpre- 
tation of observations. We present an observed relation between the large frequency 
separation, Az/, and the frequency of maximum power, t'max- We find that Ai/ oc v'^g^, 
allowing prediction of l^v to about 15 per cent given t'max- Our result is further sup- 
ported by established scaling relations for lS.v and i^max and by extended stellar model 
calculations, which confirm that lS.v can be estimated using this relation for basically 
any star showing solar-like oscillations in the investigated range (0.5 < AZ/Mq < 4.0). 

Key words: stars: fundamental parameters — stars: oscillations — stars: interiors. 



1 INTRODUCTION 

Simple scaling relations for asteroseismic quantities have 
proven very useful when analysing stellar oscillations. 
In particular, rela tions fo r scaling th e quantities z^max 
iBrown et al.lll99ll 'l and Ai/ l|Ulrichlll986l ') from the Sun have 
been widely used to constrain stellar global parameters, or 
simply to verify the asteroseismic signal of stars that al- 
ready had well-constrained stellar parameters. Here, t'max 
is the frequency of maximum power of the oscillations, and 
Ai/ is the sol-called large frequency separation between con- 
secutive overtones. For the Sun, i^max — 3100 /iHz and Ai/ 
~ 135 /xHz. 

Over the past decade, detections of solar-like oscil- 
lations in o ther stars have b een growin g in number (se e 
reviews by lAerts et al. I l2008l : [BcddinE fc KieldsenI 120081 ) . 
and the CoRoT space mission has caus ed a recent surge 
l|Michel et al.ll2008l : lde"Ridder et al.ll2009l ). This flow of data 
is ex pected to increase to a flood wit h the Kepler mis- 
sion IChristensen-Dalsgaard et al.l [ioosi ). For many stars, 
knowing Ai/ or i/m ax can help constrain the stellar param- 

[2OO8I: iKaUinger et al.ll2008l: 



tions to support this investigation, including a comparison 
of stellar models with the established scaling relations for 
Ai/ and fmax. 



2 THE OBSERVED RELATION 

Table [T] shows published measurements of A// and //max for 
55 stars including the Sun. In some cases, where the value 
of i/max was uot Specified, we have estimated it from the 
published power spectra. The location of the stars in the 
H~R diagram are shown in Fig. [T] 

In Fig. [2] we plot Ai/ versus i/max for the data listed 
in Table [T] We see a remarkably tight relation over nearly 
three orders of magnitude. A power-law fit gives 

n0.772±0.005 



Ai/ = (0.263 ± 0.009)^iHz (i/maxZ/iKz) 



(1) 



eters considerably (IStello et al 
ISteUo et all boodl : iMigl: 



lio et al 



20091 ) and hence yield im- 



portant new insights into stellar structure and evolution, 
even without a full asteroseismic analysis. In addition, re- 
cent developments of automated software for asteroseismic 
data analysis benefi t from knowing a s imple relation be- 
tween Ai/ and i^max (|Mathur et al.ll2009l : [H"uber et al.ll2009l : 
iHekker et alll2009bh . 

In this Letter we investigate the relation between Ai/ 
and i/max, building on previous studies to show how these 
two quantities are related to each other and to the global 
stellar parameters. We also carry out stellar model calcula- 



The two stars that deviate the most from this relation are 
T) Ser and ^ Hya. In both cases, there are reported ambi- 
guities in the determination of the large separation, which 
could explain the deviation. Apart from those two cases, all 
stars in Fig. [2] fall within -1-10 per cent and —15 per cent 
of the fitted relation, while the main sequence sta r s alone 
faU within ±5 per cent. We note that iHekker et all l|2009bl ) 
have confirmed the tight correlation for a sample of several 
hundred red giants observed by CoRoT, which also seems to 
indicate a tendency for some stars to fall below the relation. 



3 SCALING RELATIONS 

Can previously published scaling relations for Av and i/max 
explain the ob served relati on between them? It is well- 
established (e.g. lUlrichlligsd ) that, to a good approximation. 
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Table 1. Published measurements of Au and fr. 



Star 



(/.Hz) 



Source 



rCet 
QfCen B 
Sun 
L Hor 
7 Pav 
a Cen A 
HD175726 

Ara 
HD181906 
HD49933 
HD181420 
/3Vir 

Her 
/3Hyi 
Procyon 
?7 Boo 

Ind 
j; Ser 
«Hya 
p Vol 
eOph 
e, Dra 
fc Oph 
HR3280 

31 CoRoT giants 



170 
161.4 
134.8 
120 
120.3 
106.2 
97 
90 
87.5 
85.9 
75 
72 
56.5 
57.5 
55 
39.9 
25.1 
7.7 
6.8 
4.9 
5.3 
4.0 
4.5 
3.2 
2-7 



4500 
4100 
3100 
2700 
2600 
2400 
2000 
2000 
1900 
1760 
1500 
1400 
1200 
1000 
1000 
750 
320 
130 
90 
51 
50 
36 
35 
25 
15-73 



Teixeira et al. (2009) 
Kjcldsen et al. (2008) 
Kicldscn ct al. (2008) 
Vau clair et al. (2008) 
Mosser et al. (2008) 
K ieldsen et al. (2008) 




Garci^_et_aL 
Bouch^_e^_aL 

Garcia et al. 
GarciaetaL 
Garcia et al. 
Carrieretal. 
Bon anno et al. (2008) 
Kie ldsen et al. (2008) 
Arc ntoft et al. (2008 ) 
Carrier_et_al^ (2005a ) 
Kieldsen et al. (2008) 
Barb an et al. (2004) 
Fran dsen et al. (2002) 
unpublished WIRE data 
Barban et al. (2007) 
unpublished WIRE data 
unpublished WIRE data 
unpublished WIRE data 
Kallinger et al. (2008) 
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Figure 1. H-R diagram showing all stars listed in Table [T] The 
ellipse indicates the approximate location of the 31 CoRoT red 
giants. The evolutionary tracks (grey curves), illustrate the range 
in mass and evolutionary state that we investigate with our stel- 
lar models (see Sect.|4]l. The dashed lin es show the approximat e 
location of the classical instability strip dSaio &: Gautschvll 1998t) . 
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Figure 2. Observed Av versus i-'max for the stars listed in Table[T] 
The 31 CoRoT red giants are plotted with plus symbols. The solid 
line is a power-law fit, and dotted lines show -|-10 per cent and 
— 15 per cent deviations. 



Av is proportional to the square root of the stellar density: 



(M/M0)°-5(r<3ff/reff,0)^ 



A^^0 Vp0 {l/Lq^-^^ ■ 

Also, foUowing lBrown et al.l (|l99ll ) and lKield sen fc Bedding 
(| 19951 ), we expect i^max to scale as the acoustic cut-off fre- 
quency, Va,c- Hence, 



M/Mq(J,s/T,«.q) 
L/Lq 



(3) 



where it is observed for the Sun that i^ac.w — 1.7i^max,Q 
l|Bahnforth fc Goughl [l990l : iFossat et al]|l992l ). Given that 
Al/ and v-max scale differently with stellar parameters, the 
tightness of the correlation in Fig. [2] is perhaps surprising. 
To understand this, we raise Eq. ([Sjl to the power a and 
divide by Eq. ([2]) which, after rearranging, gives 



"(T'cff/T'cff,©)^ 



(L/Le)o.' 



^max , 



(4) 



The tight relation in Fig. [2] leads us to conclude that the 
scaling factor in square brackets in Eq. (|4]) must be approx- 
imately constant when a — 0.77. Indeed, this is easy to see 
because that value of a eliminates almost completely the 
parameter that varies the most over the stars being consid- 
ered, namely L, and leaves a very weak dependence on M 
and Tcff . This is illustrated in Fig. |3l which shows the scal- 
ing factor for the observed value a = 0.77. The scaling factor 
indeed covers a band 25 per cent wide for the region sam- 
pled by the stars (see Fig.[T|, in agreement with the -1-10 to 
— 15 per cent deviation shown by the observational data in 
Fig. [2] Figure [3] further supports the lower scatter observed 
for the main sequence stars and the higher, slightly skewed, 
scatter for red giants. We have determined the value of a 
that minimizes the spread in the scaling factor and found 
that it depends slightly on the stellar parameter range that 
we consider. For example, consideration of the entire model 
grid requires a = 0.78, while considering only ZAMS models 
gives a = 0.75. 
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Figure 3. The scaling factor in square braclcets in Eq. ^ for 
a = 0.77, for the same evolutionary tracks as in Fig. [T] The hor- 
izontal dotted curve is the ZAMS and the dashed lines indicate 
the instability strip. Masses are indicated in solar units. 



4 MODELS 

We now expand our investigation of the /S.v-Vraax relation 
using two sets of stellar models. One is a dense grid, compris- 
ing over a million mod els derived u sing the Aarhus stellar 
evolution code ASTEC (jChristensen-Dalsgaard 2008b|) and 
the ad iabatic pulsation code adipls ( Christensen-DalsgaardI 
l2008al ). This grid is slightly expanded in parameter 
space but otherw ise identical to the one constructed by 
IStello et al.l (|2009l ) (see references he rein) . The other is a set 
derived using the Yale code yrec (|Demaraue et al] |2008| ) 
with mixing-length parameter a = 1.80, initial hydrogen 
abundance X = 0.707, and metallicity Z = 0.018. Although 
the set of yrec models is more sparse, with a total of 92 
models, it spans almost uniformly the same range in pa- 
rameter space as the ASTEC grid. The region covered by the 
models is indicated by the evolutionary tracks in Fig. [1] 

In addition to the global stellar parameters M, L, and 
Tcff, we calculated for each model the large frequency sep- 
aration as the inverse of the sound travel time through the 
star: 



Af = 



(5) 



where c is the sound speed (|Tassoulll 19801 : ICoughll 19861) . We 
also calculated the acoustic cutoff fre quency by assuming 
an iso thermal atmosphere, which gives (|Balmforth fc Goughl 
Il990l ) 



c 
2H' 



(6) 



evaluated at the surface (T — Toff). Here, H = p/{gp) is 
the density scale height, p is pressure, g is gravity, and p is 
density. 

In Fig.Uwe show a subset of the ASTEC grid for a fixed 
metallicity {Z = 0.014). Contours of constant [2 J dr/c]~^ 
(magenta) are almost parallel to contours of constant c/2H 
(cyan), indicating a strong correlation. Hence, knowing one 
of these two quantities gives a good estimate of the other. 
This is particularly pronounced in the lower-right corner of 
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Figure 4. H-R diagram of grid subset (astec models with 
Z = 0.014). Models within a narrow range of fixed values 
of [2 J dr/c\~^ (magenta) and c/2H (cyan) are indicated. The 
dashed lines indicate the instability strip. Note that the low-mass 
models have been evolved beyond the age of the universe. 
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Figure 5. l\u versus v^c relation for the ASTEC models. Blue 
symbols are models hotter than the red edge of the instability 
strip, while red and grey symbols are cooler. Red symbols are low 
mass (M < 1.2Mq). Solid line is a power-law fit and dotted lines 
show +10 per cent and —15 per cent deviations. 



the diagram, corresponding to cool main-sequence and sub- 
giant stars. A somewhat weaker correlation is seen in the red 
giant phase. Note that we neglected the post He-core burn- 
ing phase, including any mass loss associated with red giant 
branch evol ution. If included, it would blur the contours in 
this region l|Stello et al.ll2008h . 

Figure [S] shows the relation between /S.u and v^c for 
the ASTEC models. We find very similar results for the yrec 
models. Stars hotter than the red edge of the instability strip 
are generally not expected to exhibit solar-like oscillations 
and are shown in blue, while cool models are shown in red 
and grey. We compare the models with observations by fit- 
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ting a power law with the observed value of a = 0.77 fixed. 
We see excellent agreement with observations, particularly 
at low masses (< 1.2 M0, red), which is perhaps not surpris- 
ing since the observations are predominantly from low-mass 
stars. Interestingly, all cool main-sequence models - high 
values of i-'ac - show a very tight relation, while for more 
evolved models - lower values of Va.c - the higher mass ones 
(grey) scatter more and fall below the power-law relation. 
This agrees with the observations and with the scaling factor 
in Eq. (|4]) (see also Fig. [Sjl . Finally, we found that the the- 
oretical IS.v-Vac relation was not notably sensitive to metal- 
licity. However, it remains to be explained why /S.v, which 
depends on the stellar mean density, correlates so strongly 
with fac, which depends on the local conditions in the at- 
mosphere. 
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5 DISCUSSION 

Several points should be kept in mind when comparing the 
model calculations with the observations. Firstly, we do not 
measure v^c directly in stars and instead rely on its relation 
to i-'max (see Eq. [3}, assuming that i^max/i/'ac does not vary 
from the solar value. In addition, i^max can be difficult to de- 
fine and measure, especially for stars with broad or double- 
humped envelopes, as seen in Procyon (|Bedding fc KieldsenI 
I2OO6I : lArentoft et al.l [2OO8I V On the other hand, Ai/ varies 
with frequency and mode degree, and its measurement from 
the power spectrum gives an average that will not be exactly 
equal to [2 J dr/c\~^. That said, Fig. [5] clearly supports the 
tight Az^-i/max relation that we observe. 

The preceding comments relate to measuring /S.v and 
J'max from observations. We should also consider the vari- 
ous ways in which these quantities can be estimated from 
the models. So far we have used Eqs. (O and Here, we 
derive the large frequency separation and the acoustic cut- 
off frequency from the models in other ways and investigate 
any significant systematic differences. The large separation 
was derived by fitting to 11 consecutive radial orders around 
fmax, which we denote Ai/flt. The result for the YREC models 
(Fig. ini left panel) shows agreement within a few per cent 
between [2 J dr/c]~^ and Ai/flt, although we note that the 
ratio [2 J dr/c]^^ /Az/gt is always higher than unity. We see 
similar results for the ASTEC models. 

We then tested two additional ways to calculate the 
acoustic cutoff frequency. Firstly, H was derived from 
the actual density gradient in the model, but still us- 
ing Eq. ((SJ. Secondly we used the full expression ;/a 
(c/2i/)\/l - 2dH/dr (Bal mforth 



the full express io 
fc Goughlll990l ). 



The re- 



sults are shown in Fig. |6] (right panel). Apart from a few hot 
models that deviate by up to 25 per cent, the agreement is 
again within a few per cent. For our purpose of investigating 
the relation between l^v and i^ac, a few per cent difference 
in either quantity is not important, but they should be kept 
in mind for other applications. 

Finally, it is interesting to investigate how well the scal- 
ing relations for /S.u and //ac (Eqs. [5] and |3J agree with 
the model calculations (Eqs. [5] and (Gjl. Figure [7] exam- 
ines how well [2 J dr/c]^^ scales with y^. For cool mod- 
els the agreement is particularly good, almost independent 
of the evolutionary state, while it deteriorates slightly for 
hot models (those below the blue line). However, we note 



Figure 6. Left panel: Ratio between [2 J dr/c\~^ and a linear fit 
to model frequencies, Ai/flt- Right panel: Ratio between C/2H2 
and c/2Hi where H2 = —dr /din p and Hi = p/{gp) (triangles), 
and between {c/2H\)sJl — 2dH/dr and c/2H\ (diamonds). Both 
panels show results for YREC models. Colour notation follows that 
of Fig. O 
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Figure 7. Ratio of Av between model calculations and solar 
scaling for ASTEC grid {Z = 0.014). Blue dashed curve shows 
models with T(,ff = 6400 K. Annotation follows that of Fig. |3] 



that [2 J dr/c]~^ systematically overestimates the large fre- 
quency separation for cool models, corresponding to 2-3 per 
cent for the Sun. 

In a similar way Fig. |8] shows that the model calcula- 
tions of c/2H follow the MT^s^'^ /L scaling from the Sun 
quite well. The ratio is always below unity because we eval- 
uate c/2H at the surface defined as T = Teff, where c/2H is 
slightly below its maximum value, which occurs further out 
in the atmosphere. As with l\v, we conclude that the best 
agreement is for the cool models. We find similar results for 
the YREC models. 



The relation between Au and i/^ax 5 




100 1000 10000 

[c/ZH] (/xHz) 



Figure 8. Ratio of i/ac between model calculations and scaling 
for ASTEC grid {Z = 0.014). Annotation follows that of Fig. [7] 



6 CONCLUSIONS 

This Letter points out that the ratio Aiy/u^ll^ varies very 
weakly with steUar parameters, so that it is essentially con- 
stant. Hence, if either one of these parameters is measured, 
this gives a very useful and robust estimate of the other 
without any prior knowledge of the stellar global parame- 
ters, L, M, Teff, and Z. We anticipate this relation can be 
used to establish the most plausible large separation in case 
of ambiguity for datasets where !^max can be determined. 
This has already b een implemented for automated analy- 
sis of Kepler data JjOekkeT et al.|[2009al : iHuber et al]|2009l : 
iMathur et al] 120091 '). In addition, we showed that the well- 
used scaling relations for Av and Vmax agree within a few 
per cent with stellar model calculations for cool models (Teff 
< 6400 K), from the main sequence to the red giant branch, 
with a slightly increasing deviation for hotter models. 
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